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"^ ' Abstract 

JL , We propose a method for evaluating the slope (and higher derivatives) of the Isgur- 

1^ ■ Wise function at the zero recoil point using lattice simulations. These derivatives are 



required for the extrapolation of the experimental data for B — > 0*1^ decays to the 
zero recoil point, from which the Vcb element of the CKM-matrix can be determined. 

The ARGUS |]I|] and CLEO [0 collaborations study the semi-leptonic decays B -^ D*li?, and 

present results for |V^b|^(co'), where ^(tu) is the Isgur-Wise function [Q, u is the dot product 

of the four- velocities of the B and D* mesons {u = vb ■ vd*), and Vcb is the cb element 

of the Cabibbo-Kobayashi-Maskawa matrix. The presentation of the results in this form is 

suggested by the "Heavy Quark Effective Theory" (HQET) (for a comprehensive review and 

references to the original literature see ref. |^). The Isgur-Wise function is normalised to 1 

at the zero recoil point, i.e. ^(1) = 1, so the extraction of the Vcb matrix element requires 

an extrapolation of the experimental results from tu > 1 to tu = 1 Q. This extrapolation 

^In principle Vcb can also be determined by fitting the data to a theoretically predicted function of uj, 
however we wish to make strong use of the rigorous result ^(1) ~ 1. 
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requires a model function or a theoretical prediction as a guide. Recently, the first lattice 
computations of ^{uj) have been performed (for discrete values of u; > 1 [^, |^, and have been 
used to extrapolate the experimental data to cj = 1, and hence to determine Vcb- In this 
letter we point out that it is also possible to compute directly the slope of the Isgur-Wise 
function at the zero-recoil point (^'(1)), as well as some higher derivatives. 

The calculations suggested in this letter complement those already proposed in ref. 0. In 
ref.0 the procedure for computing the slope directly in the HQET was explained. Here we 
discuss the computation for QCD, i.e. for heavy quarks with a finite mass. The advantage 
of using quarks with a finite mass is that there is no uncertainty in the normalisation of 
the lattice operators used to compute the slope and higher derivatives, in contrast with 
the situation in the HQET |^ . Moreover the derivatives are computed in terms of ratios of 
correlation functions, and it is likely that there will be some cancellation of systematic errors 
and statistical fiuctuations in these ratios. 

In order to determine the Isgur-Wise function it is convenient to evaluate the elastic matrix 
element {P{p')\Q'~f^Q\P{p)) , where P{p) represents a pseudoscalar meson of mass m, con- 
sisting of a heavy quark Q and a light antiquark. This matrix element can be parametrised 
by a single form-factor /(g^) defined by 

M^ = {p{p')\QYQ\P{p)) = {p + pTfiq') (1) 

where q = p — p' and q^ = 2m? {1 — uj). The Isgur-Wise function, ^{oo) is given by 

/ {2m\l - uj)) = (1 + f^{u) ) i{u) = e» (2) 

where /3(u;) represents radiative corrections which will be discussed later (/5(1) = for the 
degenerate case we are considering here), and we neglect corrections of Oll/m). 

The derivatives of the Isgur-Wise function at the zero recoil point can be expressed in terms 
of the derivatives of M^ w.r.t. the momentum p as will be explained below (in the following 
we will take fi = 4, which is a convenient choice). For simplicity consider the case where 



p' = (m, 0) and p = {po, 0, 0,^3), with po = Jm'^ + pi B ^^ ^^^^ ^^ i^ ^ function of p|. The 
derivatives of M'^(p|) w.r.t. p\ are as follows: 

MMO) = ^l,,.„ = -(«!) + ^) (3) 

M'"(0) = I^I.M = i ff (1) 4) (4) 



d(plf">'-" 2m? V ^ ' 2 



^In actual calculations we envisage averaging over results obtained with non-zero components of p in each 
of the 3 directions. 



and similarly for the higher derivatives 0, where the "s on ^ denote derivatives w.r.t. uj. 

M^ (0) and M^ (0) can be obtained from suitable combinations of two- and three-point 
correlation functions and their derivatives. Let us consider: 

C2(t.;i?3) = Ee*'""'<0|'^p(^)4(0)|0> (5) 

.g-Po*- + ... (6) 
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and 



x,w 

= e-Po'^e-^^^-'-^^^M^fps) + ••• (8) 

Ampo 

where Jp (Jp) is an interpolating operator which can annihilate (create) the heavy-light pseu- 
doscalar meson P, and V^ represents the fourth component of the vector current ZyQj^Q Q. 
Zy is the renormalisation constant relating the lattice vector current to the physical one. 
Z = |(0| Jp(0)|P)| and the ellipses in eqs.(|^) and (H) represent terms which are exponentially 
suppressed at large time separations, i.e. when t^ and t^ — t^ are both large. We assume 
that we can study the correlation functions at sufficiently large separations to neglect these 
terms. 

In eqs.(|)-(|) we have implicitly assumed that the interpolating operators Jp are local, so 
that the corresponding matrix elements Z are independent of ps. In many simulations this 
is not the case, extended interpolating operators are used to increase the overlap with the 
ground state and the corresponding Z's depend on momentum. For the sake of clarity we 
start our discussion assuming local interpolating operators, and generalise the results to 
extended operators towards the end of this letter. 

We expand the correlation functions as power series in p^: 

C2{t,;P3) = ci'\Q + ci'\Qpl + ci'\Q{plf + --- (9) 

C3(t.,t»;P3) = C^'\tx,t^)+cP{U,t^)pl + d3''\tx,t^){piy + --- (10) 

■^Throughout this letter we give the exphcit formulae for the first two derivatives. 

''Of course one could also use the lattice "conserved" vector current, i.e. the current which is conserved 
for quarks which are degenerate in mass. 



where 



ci'\Q = -^E^3(o|Jp(x)4(o)|o) (11) 

X 

Ci'\t.) = ^E^3(0|Jp(x)4(0)|0) (12) 

X 

Cf\t^,t^) = -W^l{^\Jp{w)V\x)Jl{m) (13) 



x,w 



Ct\tx,t^) = ^E^3(0|Jp(^)V^'(x)4(0)|0) (14) 
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On a periodic lattice we propose taking x| = min{x|, (L3 — 0:3)^}, where L3 is the length 
of the lattice in the z-direction, and similarly Xg = min{a;3, (L3 — Xs)^}. The derivatives 
of the Isgur-Wise function at the zero recoil point are given simply in terms of the ratios 
i?2 3 = C'23 /C'23 (where we have suppressed the time labels tx and tw). 



1 ^ ^(2) d(2) 



' 'r(l)-i) = Rf + {R?'r-R?-R?'R? (16) 
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Since ^'(1) and ^"(1) are obtained from ratios of correlation functions, it is not necessary to 
know the normalisation constant Zy 0. 

In order to obtain the derivatives of the renormalisation group invariant Isgur-Wise function 

^{uj) (defined in detail in ref. [^]) from the results for C,' and ^", it is necessary to evaluate the 

derivatives of i3{uj), the radiative corrections in eq.(0). This is a calculation in continuum 

perturbation theory, although the results depend on the mass of the heavy quark being used 

in the simulations. The ingredients necessary to perform the evaluation of /?'(1) and /?"(1) 

for elastic matrix elements are presented in the appendix of ref.P]. In the recent simulation 

by the UKQCD collaboration rriQ was about 1.36 GeV ^, for which, taking Ajg^ = 0.25 

GeV, we find /3'(1) — —0.24 and /3"(1) — 0.17 (for Uf = 3, where Uf is the number of quark 

flavours in the evolution equations, however the dependence on the Uf is mild, e.g. the results 

are -0.23 and 0.17 for Uf = 0). Having obtained the derivatives of ^{uj) at the zero recoil 

point, it is possible to obtain those of the form-factors for B ^ D and D* semi-leptonic 

decays by using continuum perturbation theory and the heavy quark symmetry]^]. 

^Zv can however be determined with excellent precision (to about 0.1% with O(IOO) gauge field config- 
urations on typical lattices) using the relation: Zy = C2{tx + ^-m; 0)/C3{tx, t^; 0). 



The discussion in this letter can readily be generalised to computations which use extended 
interpolating operators. In this case the matrix elements Z become functions of p|. Expand- 
ing ^(pg) as a power series in pg, 

Z(p^) = Z(°) + Z^^^vl + Z^'^vl + ■ ■ • (17) 

we define w^ = Z^"^'^' /Z^^\ Equations (pIS]) and (|1^) become modified to 

' 'm + D = Rf'-R?' + Rf (18) 
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' 'f'(i)-i') = /2f + (4^))^-4')-i?ri?r + 4n 
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B^^Rf - R?^Rf (19) 

Thus in order to obtain the derivatives of the Isgur-Wise function using extended interpo- 
lating operators 
as functions of t^. 



fating operators we need to determine R^ and Rz . This is done by fitting i?2 and R2 



and 



^'' = A« - (-fif ' - A) '- + 24*' + («?')' - ^^' + A (21) 



In this letter we have proposed a method for the direct determination of the slope and 
higher derivatives of the Isgur-Wise function at the zero recoil point. We expect that this 
technique, particularly when combined with standard lattice determinations of C,{uj), for 
discrete values of cj > 1, will provide a valuable tool for the extrapolation of experimental 
results to u; = 1, and hence will lead to a more precise determination of Vet- The method 
can also be generalised to other decays, such as the semileptonic decays of heavy baryons. 
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